Doppler Modulation of X-ray lines in Cygnus X-3 
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ABSTRACT 

We measured Doppler shifts of three bright spectral lines in the X-ray emission from Cygnus 
X-3 as recorded by the Chandra X-Ray Observatory. Doppler shifts of lines associated with 
Si XIV and S XVI exhibit orbital modulation. The magnitude and phasing of this modulation 
relative to the orbital ephemeris indicate the location of the source of this emission within the 
wind emanating from the compact object's companion. These observations enable us to make 
an indirect measurement of the separation of the two stars. Under certain assumptions our 
observation of a line associated with Fe XXV also limits the mass of the compact object Mq < 
3.6M . 

Subject headings: X-rays: binaries — binaries:close — stars: individual (Cygnus X-3) 



1. Introduction 

In spite of near continuous study at many wave- 
lengths since its discovery (Giacconi et al. 1967), 
the nature of the Cygnus X-3 binary system re- 
mains mysterious. Several important questions re- 
main about the nature of the compact object and 
its companion and about the specific location of 
the source of observed X-ray and infrared emis- 
sion. Existing models describe the system either 
as a high-mass system consisting of a black hole 
and a Wolf-Rayet star (Cherepashchuk & Mof- 
fat 1994) or as a low-mass system consisting of 
a neutron star and a degenerate companion (Ta- 
vani, Ruderman, & Shaham 1989). Both of these 
models account for the intense X-ray emission, the 
existence of a Helium-rich stellar wind, and the re- 
sulting large mass loss from the system. 

The observation of a 12.6 ms pulsar signal in 
gamma ray data (Chadwick et al. 1985) is evi- 
dence that the compact object is a neutron star. 
However, though this pulsar signal has been re- 
peatedly observed by the group that first reported 
it (Bowden et al. 1992), it has not been confirmed 
by independent investigators and it has never been 
observed at lower energies. Non-observation at 
lower energies is not incompatible with observa- 
tion at higher energies since the scattering dis- 



tance for X-rays in the dense wind is much shorter 
than for gamma rays. If true, the pulsar detec- 
tion is proof that the compact object is a neutron 
star. Absent independent confirmation, however, 
the reports of a pulsar signal must be regarded 
with skepticism. 

One key to understanding the nature of the two 
stars in the system would be a direct measurement 
of their masses. By measuring the Doppler shift 
of the infrared emission, Schmutz et al. (1996) re- 
port that they have measured the velocity of the 
infrared companion and constrained the compact 
object to be, at the least, a very massive neu- 
tron star. Hanson, Still, & Fender (2000) use 
the Doppler shift of a different infrared feature 
to make a more correct measurement of the ve- 
locity of the companion which prefers lower mass 
systems except that it also implies that the com- 
panion must be the dominant mass in the system. 
This rules out systems consisting of a neutron star 
with a white dwarf companion. Knowing the ve- 
locity of only one of the stars of a binary system 
does not constrain the total mass of the system. 
To determine the masses of both stars, we must 
make a measurement of the velocity of the com- 
pact object. 
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2. Observations 

We used data collected by the Chandra X-Ray 
Observatory [CXO) during the radio quiet period 
prior to the April 2000 radio outburst (McCol- 
lough ct al. 2000). Data were taken on two occa- 
sions two days apart and each observation lasted 
approximately five hours, completely covering one 
orbital period of the source. The observations 
were carried out using the High Energy Transmis- 
sion Grating Spectrometer (HETGS). Due to pho- 
ton pileup in the Advanced CCD Imaging Spec- 
trometer caused by the brightness of Cygnus X- 

3, the automatic processing was unable to pin- 
point the center of the zeroth order maximum of 
the grating spectra. We refined the automatically 
determined source location by fitting prominent 
spectral features in both the positive and negative 
first order spectra of the Medium Energy Grat- 
ing (MEG) and the High Energy Grating (HEG). 
The source location was adjusted so that the fit- 
ted wavelengths for all four first order spectra were 
identical. By this method, the location of the X- 
ray source was determined to be 20 h 32 m 25!8±0!l 
+40°57'28'.'0±0'.'l which, given the Of.' 6 absolute 
astrometric accuracy of the CXO , is consistent 
with the radio location given by Ogley, Bell Bur- 
ncll, & Fender (2001). 

We divided the data into eight bins by the 
orbital phase of the source according to the 
ephemcris in Table 1. In each of these phase 
bins, we fit a continuum model and Gaussian line 
profiles to data from both positive and negative 
first-order spectra from the HEG and the MEG 
simultaneously. The continuum was fit by a power 
law spectrum with a model of interstellar absorp- 
tion. Though the distance to Cygnus X-3 is still 
poorly determined, the degree of the absorption fit 
by the spectrum, 8.5±0.1xl0 22 H I/cm 2 is con- 
sistent with other measurements of the distance 
to Cygnus X-3 or of the intervening absorption 
column (Scrlcmitsos et al. 1975; Dickey 1983; 
Prcdchl ct al. 2000). Our measured value is also 
independent of orbital phase. 

We fit a Gaussian line function to the three 
most prominent lines visible in the spectrum, the 
Lyman a lines of Si XIV and S XVI and the He- 
lium a line of Fe XXV. Only these three lines were 
strong enough in the phase-dependent spectra that 
they could be fit reliably for all phase intervals. 



Examples of these fits are shown in Figure 1. The 
Si XIV and S XVI lines showed a consistent aver- 
age Doppler shift corresponding to a recession ve- 
locity of ~ 550 km/s. This is somewhat less than 
the velocity of 750-800 km/s reported by Paerels 
ct al. (2000), also from HETGS data, but nearly 
consistent given the absolute spectral accuracy of 
the HETGS is on the order of 100 km/s. 

Figure 2(a) shows the Doppler shift of the 
6.18 A Lyman a line of Si XIV and Figure 2(b) 
shows the Doppler shift of the 4.73 A Lyman a 
line of S XVI. The parameters of the best fit of a 
cosine function to the Doppler shift of each line are 
given in Table 2. The sinusoidal function provides 
a better fit to the data than does a constant func- 
tion. For Si XIV, x 2 P er degree of freedom drops 
from 2.7 to 1.6 with the change from a constant 
function to a sinusoidal function. In S XVI the 
change in \ 2 per degree of freedom is from 2.8 to 
1.5. The F-test probabilities for these detections 
are 18% and 14%, respectively, so the chance prob- 
ability of the simultaneous detection of these two 
modulations is 2.5%. 

The Doppler shift of the Si XIV line has signif- 
icant day-to-day variation though the modulation 
is clearly present on both days. The modulation of 
the S XVI line appears to be particularly repeat- 
able. Sinusoidal functions may not be the most 
appropriate functional forms to fit to these data. 
In the simplest configuration, the emission region 
will be in the form of a ring or a disk near the 
compact object sampling different wind velocities. 
In addition, different parts of the emission region 
may be occulted at different times. Beyond the 
apparent inadequacy of the shape of the fit func- 
tion, the phase of the maximum Doppler shift may 
be affected by the differential occultation of the 
emission region. Nevertheless, the amplitude of 
the fit sinusoid gives a reasonable approximation 
of the amplitude of the correct functional form. 
The quality of these fits to the data (x 2 per degree 
of freedom is 1.6 for Si XIV and 1.5 for S XVI) in- 
dicate that within our uncertainties the true func- 
tional form is not much more complicated than a 
sinusoid. 

Figure 2(c) shows the Doppler shift of the 
1.86 A Fe XXV He a line. The data are consis- 
tent with a constant Doppler shift. The x 2 per 
degree of freedom assuming a constant Doppler 
shift is 1.77 while the best fit sinusoidal function 
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Table 1 

Cygnus X-3 orbital ephemeris 



T n = T Q + P n + P Pn 2 /2 


To 


=JD 2440949.89016 ±0.00064 


Po 


=0.19968462 ±6 x 10~ 8 d 


P 


= (5.52 ±0.12) x 10~ 10 



Note. — T n are the times of suc- 
cessive X-ray minima from Stark & 
Saia, in preparation. 



Table 2 

Doppler Shift of Strong Emission Lines 



Line 


Source 
Wavelength 

(A) 


Observed 
Wavelength 

(A) 


Recession 
Velocity a 
(km/s) 


Modulation 
Amplitude 

(A) 


Modulation 
Amplitude 

(km/s) 


Phase 
Offset b 

(°) 


Si XIV 
S XVI 
Fe XXV 


6.1822 
4.7292 
1.8617 c 


6.19381 ±0.00058 
4.73774 ±0.00045 
1.86340 ±0.00038 


563 ±28 
541 ± 29 
274 ± 61 c 


0.0027 ±0.0007 
0.0025 ± 0.0006 
< 0.0014 d 


133 ± 35 
158 ±35 
< 220 d 


21 ±14 
-37 ± 14 



a Quoted uncertainty is from the fit of a Gaussian line profile. Absolute spectral accuracy of the HETGS 
is on the order of 100 km s . 

b Phase offset is for the maximum red-shift relative to the time of the minimum intensity of the X-ray 
continuum. 

The reported Fe XXV source wavelength is from a combination of lines we may have modeled incor- 
rectly. The value of the recession velocity is, therefore, subject to a large, unknown systematic error. It 
is not, however, important to the results presented here. 

d Upper limits are 90% confidence limits. All other uncertainties are la 
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Fig. 1. — Examples of line fitting from this anal- 
ysis. The fits indicated here by the thick line are 
for the phase interval centered on 0.625 on 2000 
April 4, The left panels are for FE XXV, the cen- 
ter panels are for S XVI and the right panels are 
for Si XIV. The top panels are from the positive 
first order HEG spectrum and the bottom panels 
are from the positive first-order MEG spectrum. 
The fitted profiles are derived from fitting all four 
of the first-order spectra simultaneously. 
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Fig. 2. — Observed wavelength of (a) the Si XIV 
Ly a line and (b) the S XVI Ly a line and (c) 
the Fe XXV He a line as a function of the X-ray 
phase. The S XVI point at phase 0.25 from 2000 
April 6 is corrupted by a peculiarity of the back- 
ground and is excluded from the fit. The phases 
on the plots have been shifted slightly to reveal 
overlapping data points. The dashed line in (c) 
corresponds to the 90% confidence upper limit on 
the modulation of the Fe XXV line. 



yields a x 2 per degree of freedom of 1.84. The 
non-observation of a modulation of the Doppler 
shift in Fe XXV allows us to set an upper limit 
on the velocity of the emission region for this line. 
The 90% confidence upper limit derived for a fit 
to these data is 220 km/s. 

3. Discussion 

In order to use our observations to draw conclu- 
sions about the size and mass of the Cygnus X-3 
system, we must first decide where in the system 
the emission is produced, van Kerkwijk (1993) 
observed that line emission from hydrogen-like he- 
lium was produced in that part of the stellar wind 
of the companion which is shaded by the compan- 
ion from the X-rays emanating from the compact 
object. According to his model, the rest of the 
wind is too highly ionized to produce helium line 
emission. The coincidence between the maximum 
red shift of the line emission from hydrogen-like 
silicon and sulfur and the time of the X-ray mini- 
mum (Taken to indicate the superior conjunction 
of the compact object) suggests that these emis- 
sion features are produced in the companion wind 
near the compact object. The fact that the sili- 
con and sulfur emission are produced much nearer 
to the X-ray source than the hydrogen emission 
suggests that the dominant source of ionization in 
the Cygnus X-3 system is the X-rays. The level of 
ionization appears to increase with proximity to 
the source of X-rays. Emission from hydrogen- 
like and helium-like iron will thus be produced 
even closer to the X-ray source — presumably the 
compact object — than is the silicon or sulfur emis- 
sion. We will, therefore, assume that the emission 
feature produced by helium-like iron is produced 
close to the compact object. 

3.1. Stellar Mass 

The emission features of highly ionized iron 
may be produced in a region of the wind cap- 
tured by the compact object or from an accretion 
disk around the object if not from the surface of 
the object itself. Without knowing the phase rela- 
tionship between the Doppler shift of the Fe XXV 
line and the X-ray ephemeris, we can not be sure 
that it would reveal the motion of the compact 
object. If, however, we assume that this is so we 
can use our non-observation of Doppler shift in the 
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Fe XXV line to constrain the velocity of the com- 
pact object. The equation for the mass function 

m/ = ^P< ax (l-e 2 ) 3/2 (1) 

yields an upper limit on the mass function of 
0.22 M Q taking P = 4.8 h and our 90% confi- 
dence upper limit of 220 km/s for the maximum 
radial velocity i>Max- Ghosh et al. (1981) show 
that a high mass-loss rate can drive the system 
into an eccentric orbit as well as cause it to circu- 
larize. For simplicity we will assume e = for the 
time being. 

By measuring the Dopplcr shift of the cen- 
troid of an infrared absorption feature in the stel- 
lar wind Hanson ct al. (2000) calculated a mass 
function for the wind producing companion of 
0.027 Mq. We can combine this mass function 
with our mass function upper limit by simultane- 
ously solving the following equations 

m (M c sinz) 3 m _ (Mpsinz) 3 

m/D (M D + M c ) 2 ' m/C (M C + M D ) 2 

Where the subscript C refers to the compact ob- 
ject and the subscript D refers to the mass do- 
nating companion, i is the inclination of the bi- 
nary orbit. The upper limits on the masses of the 
two stars derived this way are Mq < O.24M and 
M D < 0.49M Q for an inclination, i = 90°. These 
limits essentially exclude all stars that could pro- 
duce the observed X-ray flux. Mass upper limits 
as a function of orbital inclinations are shown in 
Figure 3(a). 

The orbital model of Cygnus X-3 that has been 
used most extensively is one of two stellar wind 
models proposed by Ghosh et al. (1981) which 
have a large orbital inclination and a large eccen- 
tricity. Even Hanson ct al. (2000) state that the 
inclination should be large, but our results, taken 
in combination with theirs, make large inclinations 
unphysical. In light of other recent work (Singh et 
al. 2002) it appears that the lack of apsidal mo- 
tion in the system requires that the orbital eccen- 
tricity be small. In addition, several authors (van 
der Klis & Bonnet-Bidaud 1989; Mioduszewski 
et al 2001) demonstrate that the orbital inclina- 
tion of the system must also be small. The best 
fit to the Ghosh et al. (1981) model based on a 
variable luminosity X-ray source inside its com- 
panion's stellar wind satisfies the need for both a 



small orbital inclination (i = 24°) and eccentricity 
(e = 0.14). Using these values in Equations 1 & 2 
gives an upper limit for the mass of the compact 
object, Mq < 3.6M and for the companion of 
Md < 7.3M . These mass limits don't preclude 
the possibility of the system being composed of 
a black hole and a Wolf-Rayet star. Neither do 
they rule out the compact object being a neutron 
star but the suggestion that its companion is more 
massive than the compact object rules out the pos- 
sibility that the companion is also degenerate (a 
white dwarf). 

3.2. Stellar Size 

Table 2 shows that the maximum red shift of 
the Si XIV and S XVI emission lines are nearly 
coincident with the minimum of the continuum 
X-ray intensity. The coincidence of the maximum 
radial velocity of the line emission and the X-ray 
minimum means that the line emitting region is 
moving most directly away from Earth at the same 
time that the compact object is at its furthest dis- 
tance from Earth. This, then, suggests that the 
line emission emanates from a region of the com- 
panion's stellar wind blowing by the compact ob- 
ject and close to it. We can therefore use our mea- 
sured velocity for this wind with a model of wind 
acceleration to constrain the size of the binary sys- 
tem. 

The standard model of stellar wind acceleration 
appropriate to Wolf-Rayet stars is 

«(r)^ 00 (l-^) /3 (3) 

where r c is the core radius of the star and (3 ex- 
presses the degree of acceleration. Alternately, the 
wind velocity can be expressed in terms of the ex- 
tent of the accelerating wind, 

v(r)= Voo ^y (4) 

(Cherepashchuk & Moffat 1994) where Too is the 
distance from the wind producing star to where 
the terminal velocity is reached. We can combine 
the wind velocities implied by the analysis of the 
Si XIV and S XVI lines to provide an estimate 
of the wind velocity in the emission region sur- 
rounding the compact object of 146 ± 50 km/s. 
Taking Voo = 1500 km/s (Schmutz ct al. 1996), 
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r M = 15i?0 (Fender, Hanson & Pooley 1999), 
and 7=1 (Antokhin & Cherepashchuk 2001), 
the distance between the emitting region and the 
wind producing companion star is 1.5 ± O.5i?0. 
This result assumes an inclination i = 90°. For an 
inclination angle i = 24°, we calculate 3.6 ± 1.2_R© 
for the distance between the donor star and the 
emission region. Since this emission comes from 
near the X-ray star, the distance calculated be- 
tween the emission region and the companion star 
is the distance between the two stars. 

The size of the Roche lobe of each star can be 
calculated directly from the mass functions. The 
radius of the Roche lobes is given by 

r ± = 0.49g 2 / 3 
a 0.6 <Z 2 / 3 + ln(l + (? 1 /3) ( > 

(Eggleton 1983) where a is the binary separation 
and q is the ratio of the stellar masses, which is 
only dependent on the mass functions and which 
in the case of the results presented here Md /Mq = 
2.0. This yields a Roche lobe radius of the com- 
panion star of i?D = 1-6 ± 0.5i? Q at an inclina- 
tion, i = 24° , which is somewhat small for a Wolf- 
Rayet star (Moffat & Marchenko 1996) especially 
given that the dominant mass-transfer process is 
the stellar wind rather than Roche Lobe overflow. 
The companion Roche lobe radius is plotted as a 
function of the inclination of the system in Figure 
3(b). 

4. Conclusion 

We have determined an upper limit on the ra- 
dial motion of the X-ray emitting star. This limit 
assumes that the Doppler shift of the emission 
feature produced by Fe XXV would reveal the 
motion of the X-ray emitting star. The relative 
phases of the observed Doppler shift of the He II, 
Si XIV, and S XVI emission features suggest that 
the degree of ionization of the circumstcllar ma- 
terial is correlated with proximity to the X-ray 
source. This is consistent with our assumption. 
Combined with several other observations, we be- 
gin to paint a consistent picture of the Cygnus X-3 
system. The picture favored by our analysis is of a 
system consisting of two relatively low mass stars, 
a compact object Mq < 3.6M and its companion 
Md < 7.3-M© which orbit each other separated by 
3.6 ± 1.2i? Q in relatively circular e = 0.14 orbit 



that has a small inclination i = 24° relative to our 
line of site. Future X-ray observations will enable 
us to test the assumptions underlying this picture 
and, if they prove correct, to determine the orbital 
parameters of the system. 
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Fig. 3. — (a) Mass upper limits of the compact ob- 
ject (solid line) and the companion (dashed line) 
and (b) the Roche lobe radius of the companion 
(dotted lines are la limits) in the Cygnus X-3 sys- 
tem as a function of inclination angle. 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.0. 
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